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Abstract 

This paper investigates simple means of analyzing the error rate performance of a general q- 
ary Galois Field network coded detect-and-forward cooperative relay network with known relay error 
statistics at the destination. Equivalent relay channels are used in obtaining an approximate error rate 
of the relay network, from which the diversity order is found. Error rate analyses using equivalent relay 
channel models are shown to be closely matched with simulation results. Using the equivalent relay 
channels, low complexity receivers are developed whose performances are close to that of the optimal 
ML receiver. 

Index Terms 

Cooperative diversity; detection and estimation; source/channel coding. 

I. Introduction 

Cooperative communication is effective in extending the coverage area and increasing the 
quality of service of relay networks by means of diversity [1]. The simplest method of trans- 
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mitting a message in a cooperative communication network is routing, where the data received 
at a node are simply forwarded to another. In network coding (NC), multiple data received at a 
node can be combined and transmitted simultaneously. Hence, NC can be used to increase the 
throughput of the system and reduce bandwidth and transmission energy consumption [0. 

Non-binary NC is a simple and efficient means of achieving improved diversity gain in wireless 
relay networks compared to binary network coding. In J3]|, a q-ary Galois field (GF(q)) NC multi- 
user cooperative wireless relay network is considered. The users act as each other's relay nodes 
and use GF(q) NC and decode-and-forward (DF) techniques. The use of non-binary network 
coding was shown to lead to higher diversity orders compared to binary network coding 0. 

Amplify-and-forward (AF) and DF are two commonly used techniques for relaying. In AF 
relaying, the signal received at the relay is transmitted after being amplified by a factor. In DF, 
the received signal is decoded then re-encoded and modulated at the relay before re-transmission. 
AF needs the redesign of many transceivers while DF is computationally complex for the relay 
node. The detect-and-forward (DetF) technique, where the received symbol is detected and then 
modulated without decoding and re-encoding, is a low-complexity alternative to DF 0~||. In this 
paper, the DetF method is considered. 

Availability of channel state information (CSI) or SNR values at the relay and/or at the 
destination improves the performance of relay networks significantly. If CSI is available at an 
AF relay of a source-relay-destination (S — R — D) link, the amplification factor may be chosen 
to compensate for the instantaneous channel gain of the S — R link H. If SNR values of the 
S — R link are available at the destination in DetF relay networks, the error propagation effects 
can be mitigated by taking the relay detection error statistics into account at the destination |Q. 
In O, the closed-form error probability of a GF(2) NC DetF relay network with known relay 
error statistics is obtained. The error performance analyses in [5J and [6l are derivable when 
GF(2) NC and BPSK modulation techniques are employed, however, the analyses are highly 
complex for higher Galois fields and modulation orders. 

In order to analyze the error rate of the relay channel, various analytically tractable approx- 
imations may be used. In 0, the max-log approximation is used in the performance analysis 
of a coded, cooperative, two-hop, DetF relay network with multiple relays. The approximation 
exhibits accurate error performance for the relay network. However, the approximation in [7 J does 
not consider network-coded relay networks. Another approximation technique is the equivalent 
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relay channel method, which approximates a two-hop relay channel with a single-hop channel 
ll8l- lfT3l . The most commonly used equivalent channel model in DF and DetF relay channels 
is the minimum equivalent relay channel [|8l- [fT0ll . In this model, the SNR of the equivalent 
single-hop channel is equal to the minimum of the S — R and R — D channels' SNR values. The 
minimum equivalent relay channel model is preferable due to its simplicity. In the equivalent 
relay channel model in [[TT|. the equivalent SNR value is a function of the error probabilities 
of the S — R and R — D channels and is shown to asymptotically approach the SNR value of 
the minimum equivalent relay channel. The equivalent relay channel model is used in [fTTI to 
develop a diversity achieving cooperative maximum ratio combining (C-MRC) receiver and in 
the diversity analysis of C-MRC. However, |fTTfl is short of the diversity order and error rate 
analysis for ML detection. The equivalent relay channel model developed for multi-hop DF relay 
networks without network coding in |fTTfl. is shown to also be valid for GF(2) NC DF two-way 
relay channels in |[T2l|. In lfT3l . an equivalent relay channel method called propagation error 
modeling (PEM) is devised for a low-complexity detector for GF(q) and complex field network 
coded (CFNC) DF relay networks. In PEM, the detection errors that occur at the relays are 
modeled as virtual noise at the destination in order to overcome the error propagation effects 
in DF systems. The PEM method, which displays close performance to the relay network, can 
be used in obtaining the diversity order of the relay network. However, this method does not 
readily provide the error rate. 

To the best of our knowledge, simple means of obtaining the diversity order and error rate 
for ML detection of general GF(q) NC DetF relay networks which use relay error statistics at 
the destination is lacking. Low-complexity alternatives to the optimal receiver for such a relay 
network are also not available. The aim of this paper is to address these problems. 

In this contribution, the error rate analysis methods of (H, [fTTI are generalized to GF(q) NC 
DetF relay networks. Two-hop relay channels are approximated using equivalent relay channel 
models, where the equivalent relay channel model in ifTTI is adapted for general GF(q) NC DetF 
relay networks and is compared to the minimum equivalent relay channel model. The equivalent 
relay models are then used in finding pairwise error probabilities (PEPs), from which the diversity 
order of ML detection and an approximation to the union bound are obtained. The optimal 
receiver's error performance is compared with the ML performances of the equivalent relay 
channel receivers. It is observed that the equivalent channel models successfully approximate 
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the relay network. 

The organization of this paper is as follows. In Section II, the considered system model 
and data transmission technique are given. In Section EH, error rate analyses of relay networks 
with error-free and error-prone S — R channels are provided; the ML rule is derived and an 
approximate union bound for error rates is obtained using equivalent relay channel models. In 
Section IV, the ML rule and union bound symbol error rate (SER) performances of the equivalent 
relay networks are obtained and compared to the relay network's ML performance. The results 
are summarized in Section V. 

II. GF(g) Network Coded Detect-and-Forward Relay Model 

As a general model, the relay network model in Fig. Q] is considered. M-ary phase shift 
keying (PSK) and time division multiplexing (TDM) techniques are used in the transmission of 
the GF(q) network coded symbols. All channels are assumed to be independent and Nakagami-m 
fading. The users act as each other's relay nodes. The DetF method is used in relaying. Imperfect 
detection at the relays is taken into account, and the detection error statistics observed at the 
relays are assumed to be available at the destination for ML detection. The users and destination 
are assumed to have SNR knowledge of the incoming data's channel. 

The data of each user terminal T n , 1 < n < N, are transmitted to the destination D, directly 
and via other users as follows. T n broadcasts its modulated data to the other terminals and D, 
where the data received at the terminals are detected. After all users have completed transmitting 
their data, K — N (K > N) users transmit network coded, modulated data to D using a pre- 
determined network coding rule, hence a total of K symbols are received at D. The received 
data at D can be expressed as 

y D = H D s + n D (1) 
where y D is the K x 1 sized vector consisting of K data received at D: 

Yd = [yiD ,sys ■ ■ ■ VND,sys VlD,nc ■ ■ ■ y(K-N)D,nc\ 

(2) 

y n D,sys is the observation at D corresponding to the n th direct transmission, which is the 
systematic symbol from T n ; yw, nc is the observation at D corresponding to the I th network 
coded data, 1 < / < K — N; Hd is the K x K sized diagonal fading coefficient matrix observed 
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at D; nn is the K x 1 sized noise vector at D and s is the K x 1 sized vector comprising the 
network coded modulated data transmitted from the terminals: 

s = tp (G T u + e R ) . (3) 

Here, </?(•) denotes the constellation mapping and G is the generator matrix represents how 
network coding is performed in the relay network. For the considered network coding operation, 
G = [I : P]nxk ' wnere me N x N identity matrix I represents the systematic data (first iV 
symbol transmissions), and the (K — N) x K parity check matrix P represents the remaining 
K — N network coded symbols transmitted. The data vector sent by the user terminals is shown 
by u = [ui ... un] ' where u n is the data of T n ; e R is the K x 1 sized coded detection error 
vector at the relays; the elements of G, u and eR are in GF(q), 1 < n < N . 

Let us define the set of all modulated GF(q) network coded codewords \- When the relay 
network in Fig.[T]has error-free S—R links, there are q N possible such codewords. Corresponding 
to each source symbol configuration is a codeword X(i), where X(i) is the i th element of the 
codeword set \'- 

X(z) = <^(G T u) (4) 

assuming the i th element of the set of all data vectors of u is transmitted. The aim of the 
destination is to decide which codeword was sent based on the observation y D and the error 
statistics of the relays, where the decision at D is denoted by X(j). 
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The relay network in Fig. Q] contains multi- source relay channels due to the use of network 
coding at the user terminals. For analytical tractability in the error rate analysis, the multi-source 
relay channels will be approximated by equivalent relay channels as follows. 

Data received from the terminals T^n ^ n, are detected and then GF(q) network coded at 
the relay terminal T n , 1 < n, h < N 

UnR = UnR + e nR (5) 

where u nR is the error- free network coded data at T n ; e nR is the network coded detection error 
at T n ; u nR ,u nR and e n R are GF(q) elements. The network coded data received from T n at D is 

+ n nD<nc (6) 

where h nD nc is the fading coefficient of the T n — D channel over which the network coded data 
is sent; s n R = ^{u nR ) is the modulated network coded data sent from T n and n n £> nc is the noise 
term at the destination of the observed network coded data. When hard-decision ML detection 
is performed at the destination, the hard-decision of the data in © is used 

hard— decision 

UnD,nc — r ZnD,nc {') 

where z n D>nc is the hard-decision of y n D,nc at D. The instantaneous end-to-end symbol error 
probability (SEP) of the relay channel is: 

PrOnD.nc ^ u nR \u nR ) = p{z nD ,nc ^ u riR \u riR = u nR + e nR )p(e nR ). (8) 

The multi-source relay channel described above can be approximated by a single-hop S — D 
channel as shown in Fig.[2l where T n corresponds to the relay node R and terminals T n represent 
the source nodes S n , n ^ n. Dropping the n's, the instantaneous SEP given in ([8]) is approximated 
by the instantaneous SEP of a single-hop channel which has an instantaneous SNR per symbol 

leg- 

Peqileq) = ^{z D ,nc ^ Ur\ur) 

q -i (9) 

= p( z d ^ ur\u r = u r + e R )p(e R ). 

Two equivalent relay channel models will be used in the error rate analysis, which will be 
discussed in Section Hill 
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Fig. 2. Equivalent relay channel of a relay channel with N sources 

III. Error Rate Analysis 

Before investigating the error rate analysis, we provide the soft- and hard-decision ML rules 
for relay networks with error-free S — R links, which will later on be used in deriving the ML 
rules for relay networks with error-prone S — R links. 

Soft-decision ML is defined as ML detection at the destination based on the soft-decision 
observation vector y D in (OQ). The soft-decision ML rule of terminal Tj's data Ui, 1 < i < N, is 
given as: 

uf L = argmaxp(y D |ui). (10) 

Ui 

Using the independence properties, the likelihood function given above is expanded 

P(yr>\ui) = ^Piynlui^^piu 1 ) 

K (ID 

= ^2Ylp(yD,i\ u ^ u X u ") 

where u 1 = u\uj is the data vector excluding Ui and yo,i is the i th element of yD- Assuming q- 
ary symbols of all N users are transmitted equally likely and independently, the joint distribution 
of u 1 is: 

pin 1 ) = - . (12) 
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Inserting (fT2l) into (fTTT) gives: 



K 



p(yDk) = g- (JV - 1) X)II^A<l«i»u i ) (13) 

The hard-decision ML rule is defined as ML detection at the destination, based on the 
hard-decision observation vector z D , where each element is the individual hard-estimate of the 
corresponding element of yr>: 

hard— decision 

Yd -> z D . (14) 
The hard-decision ML rule of terminal Tj's data itj, 1 < « < iV is 

uf 41 = argmaxp(z D |«j) (15) 
where the likelihood function given in (fl"5T) is expanded using the independence properties 

K 

p(zdK) = y^TTM^K) u ')p( u1 )- (16) 



A. Relay Networks with Error-Free Source-Relay Channels 

The error rate analysis of the relay network with error-free S — R links is provided as a 
preliminary to the error rate analysis of the relay network with error-prone S — R links. 

In the error rate analysis of the relay network with error-free S — R links, the reliability of the 
observation at the destination depends only on the R — D channel SNR values. The transmitted 
block in Euclidean space is the modulated codeword in ©. The symbol error probability (SEP), 
is upper bounded by the average of the average pairwise error probabilities (PEPs) over a single 
symbol error event 031 

^<^E E p (x(») a7) 

*=i x(i),*Cj-)«xj*?«* 

where X(i) = cp(ii . . .%k) is the transmitted modulated codeword and X(j) = ip(ji . . .jk) is 
the decision at D. The codewords X(i) and X(j) are found using © and i k and j k are GF(q) 
elements, 1 < k < K. Note that each PEP P \ — > X(j) J in (fTTI) is obtained after averaging 
over the channel Hd- This averaging can be performed using the moment generating function 
(MGF) method Ifl4l [fT51: 

p (x( ^x W )^/n^(-M|i_^) d , 
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Assuming a Nakagami-m fading environment, M k (s) = (1 — s^k/m)~ m and % are the MGF of 
the instantaneous SNR and average SNR value of the k th element of the received data vector in 
(QQ), respectively and m is the fading figure. 

Obtaining the PEPs is not as straightforward, however, when the S — R channels are error- 
prone. The equivalent relay models described in the next section can be used to overcome this 
difficulty. 

B. Relay Networks with Error-Prone Source-Relay Channels 

In relay networks with error-prone S — R channels, when no relay error statistics are available 
at the destination, the soft-decision ML rule is equivalent to that used in the error-free case, 
given by (fTTI) . However, if the relay error statistics are known at D, the likelihood function for 
the soft-decision ML rule is 

p(ydK) = ^p(yDk,eR>(e R ) 

(19) 

where the fact that relay events are independent from the source symbols is used. Using (fT2~l) : 
p(YdK) = q~ {N ~ 1] ^^p(y D |u i ,u 1 ,e R )p(e R ). ( 2 0) 

eR u' 

Similarly, for the hard-decision ML rule, when the relay error statistics are available at D, the 
likelihood function is: 

p(zdK) = <T (Ar ~ 1} ^p(z D |Mi, u\ e R )p(e R ). (21) 

eR u' 

Compared to (fT3l) and (fl6l) . the likelihood functions in (1201) and (|2T|) involve the averaging over 
p(e R ), which increases the mathematical complexity of the soft- and hard-ML rules, respectively, 
hence increasing the computational complexity of the optimal soft- and hard-receivers. The 
computational complexity of the receivers can be decreased using equivalent relay channel 
models, which will be discussed later on. 

For the relay networks with error-prone S — R channels, the error rate analysis is not as simple 
as in the error- free case. This is due to two facts. First, the optimal ML receiver is complex 
because of the averaging of S — R link errors. Second, even if a suboptimal receiver is used 
which does not perform this averaging, the analysis is still complex due to the S — R error 
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events. In order to obtain a tractable analytical results for the PEPs, we propose using equivalent 
relay channel models both for the ML receiver and for its analysis. Thus, we will approximate 
the error-prone relay channel with an equivalent error-free relay channel model, and analyze the 
performance of the optimal receiver receiver for that model. 

In ifTTI . an equivalent relay channel model is proposed that approximates a multi-hop, coop- 
erative DF relay channel with a single-hop channel. The model is shown to be valid for any 
modulation order, assuming coherent modulation/demodulation. The equivalent relay channel 
takes the channel model of the S — R and R — D channels, which is exemplified for a flat 
Rayleigh fading, two-hop relay channel with coherent BPSK modulation/demodulation, as shown 
in Fig. [3l The instantaneous bit-error probability (BEP) of the equivalent relay channel is found 
as ED 

Pq-* = pb SR (1 - Pro) + p rd (l - P b SR ) (22) 

where Pg R and P^ D denote the BEP's of the S — R and R — D links, respectively. Assuming 
coherent BPSK modulation/demodulation, the equivalent BEP expressed above is approximated 
to a single-hop channel that has an instantaneous SNR per bit 'Jq-i ifTTfl: 

Pq-^q{\Jh^) (i-^(\/2 : ^))+q(v /2 7^) (^(a/ 2 ^)) 

(23) 

where r yg R and "y RD denote the instantaneous SNR's per bit of the S — R and R — D links, 
respectively. From (1231) . an instantaneous SNR value of the equivalent relay channel is found as 

CD 

7Q-> = \{Q-\P b Q-i)Y ■ (24) 

Due to the use of the Gaussian Q-inverse function in (T2"4l) . the equivalent relay channel model is 
referred to as the Q-inverse equivalent relay channel in the rest of this paper. Notice from (|22|) 
and (|23l) that, for a fading channel, the instantaneous SNRs 7^ and "f RD are random variables, 
thus the instantaneous SNR Jq-i is also a random variable. For the error rate analysis, the 
distribution of 7^_i is approximated by a conventional fading SNR distribution whose parameters 
are determined. 
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Fig. 3. Equivalent relay channel of a single-sourced relay channel 



The Q-inverse equivalent relay channel model, originally developed for DF relay networks, is 
used in a GF(2) NC DF two-way relay channel in [[T2| and has the same BEP expression as in 
(|23]l and SNR value in (|24j). 

The novel part of this paper involves the adaptation of the Q-inverse equivalent relay channel 
model to GF(q) NC, coherent M-PSK modulated, DetF, Nakagami-m fading relay channels with 
multiple sources. The proposed equivalent relay channel model is used in obtaining approximate 
error rate expressions for GF(q) NC, coherent M-PSK modulated, DetF, Nakagami-m fading 
relay networks. Soft- and hard-decision ML receivers employing the Q-inverse equivalent relay 
channel model, which display very close performance to the optimal receiver, are developed. 
The instantaneous SEP in © is approximated to the instantaneous SEP of a single channel 



(25) 



P*_ x =2Q( v /2^7sin(7r/M)) 
= Pr(z Djnc ^ u R \u R ) 

for coherent M-ary PSK for large values of M and SNR [1J. The instantaneous SNR per symbol 
of the Q-inverse equivalent relay channel is found using (1251) as: 

(Q-^O.SP^)} 2 



7q- 



2 sin 2 (vr/M) 



2 sin 2 (vr/M) 
The average value of (|26l) 



Q 1 ^0.5 ^2 P( z D,nc 7^ u R \u R = u R + e R )p(e R )j 



(26) 



2 sin 2 (vr/M) 



E 



{Q-^O.bP^)}' 



(27) 



has no known closed-form expression but can easily be obtained numerically offline and used 
in online calculations with a look-up table. 
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Fig. 4. Histogram of the Q-inverse equivalent channel's instantaneous SNR value in a Nakagami-m fading environment (m=2) 



In a Nakagami-m fading environment, it is seen in the histogram result in Fig. |4] that the 
distribution of Jq-i is well approximated by a Gamma distribution when m = 2, so the equivalent 
relay channel can be approximated as a single-hop Nakagami-m channel. The histogram is 
obtained via Monte Carlo simulation of the instantaneous SNR in (|26l ). Based on this result, the 
instantaneous SNR value in (|26|) is assumed to be Gamma distributed: 

Q r(m) V7Q-i/ V 70-i / 

The Q-inverse equivalent hard-receiver of the relay network in Fig. \T\ obtains the ML estimates 
of Ui with the hard-decision ML rule in (fl"5l) . where the SEP terms of network coded data in the 
likelihood function in (|2~D) are approximated by (T25T) . The Q-inverse equivalent hard-receiver is 
illustrated with an example as follows. Consider the QPSK modulated, 2-user GF(4) NC DetF 
relay network in Fig. Q] for iV = 2, which has a generator matrix: 

/ 1 1 1 \ 
Gi = . (29) 

V 1 2 1 / 

Hence, network coding is performed at Ti and T 2 as u\r — u\ + 2u 2 and u 2 r — U\ + u 2 , 
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respectively, where Ux, u 2 , um, u 2R G GF(4). The received data vector at D is detected as 

P ..rp hard — de cision . -.rp 

[yiD,sys V2D,sys VlD,nc V2D,nc\ ~ [Z\D,sys Z>2D,sys Z\D,nc Z 2 D,nc\ 

where y n D, S ys an d UnD,nc denote the systematic and network coded data received from T n at 
D, respectively; z n D, S ys and z nD ^ nc are the hard-decision data obtained by individually detecting 
UnD^sys and VnD^nc respectively, n = 1,2. The likelihood functions of the detected data in (1301) 
are used in the hard-decision ML rule in (TT3T) . where the likelihood functions of the network 
coded data are approximated using the Q-inverse equivalent channel method as follows. The 
instantaneous SEP of the relay channel T 2 — T\ — D, over which U\ R is sent, is approximated 
using ((25]): 

» (3D 

= 2^ P \Z\n,nc ^ Ui R \u 1R = u 1R + e 1R ) p{e 1R ) 

eiK=l 

where e 1R is the network coded coherent QPSK demodulation error at T\ and Jq-i i is the 
instantaneous SNR per symbol of the equivalent relay channel of T 2 — T x — D, which is found 

as 

7Q-Sx = W- 1 (0.5^-i, 1 )} a 

/ 3 \} 2 (32) 

Q' 1 I 0.5 ^ p (^id,t»c 7^ wi-rImia = M ifl + eifl) p( e ifl) 

Similarly, the instantaneous SEP of the relay channel Ti — T 2 — D, over which u 2 # is sent, is 
approximated 

^ (33) 

= 2^ p (^2D,nc 7^ u 2R \u2R = u 2R + e 2R ) p(e 2R ) 

where e 2R is the network coded coherent QPSK demodulation error at T 2 and 7q-i, 2 is the 
instantaneous SNR per symbol of the equivalent relay channel of Ti — T 2 — D: 

7 Q -2 = {g- 1 (0.5P^ li2 )} 2 

- / 3 \y (34) 

Q 1 I 0.5 22 p{z 2D>nc ^ u 2R \u 2R = u 2R + e 2R )p(e 2R ) \ > . 

V e 2R =l J ) 

Another equivalent relay channel model is the minimum equivalent relay channel model, 
which is commonly used in analyzing the error performance of DetF relay networks ll8l- [fT0ll . 
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The minimum equivalent relay channel has a lower computational complexity than the Q-inverse 
equivalent relay channel, which is preferable in the receiver design. In this model, the two-hop 
relay channel shown in Fig. [3] is approximated to a single-hop channel, where the equivalent 
instantaneous SNR per symbol, 7 m j n , is chosen as the minimum of the S — R and R—D channels 
instantaneous SNR values [|8l-[fT0ll 

j min = mm{-f S R,lRD}- (35) 

Here, jsr and 'Jrd are the instantaneous SNR's per symbol of the S — R and R—D channels, 
respectively. In a Nakagami-m fading environment, the SNR value in (|35l) is distributed as lfT6l 

man , \ rn RD -l / N fc 



flmin (7) 



1 (!^iY SR e -{m + ^h y 

r(ms^) 7 V Isr J trt k\ \ 7-RD / 



k ~° (36) 

-1 / \ raRD , „, \ m SR — 1 -, / \ k 

+ 1 ( n^£J] e -(^ + ^> y 1 ( 

F(m RD )i V Ird J ~ k\ V 75iJ / 

where m Sfi and m RD are the fading figures of the S — R and R — D links respectively; 7s.r 
and 7rd are the average SNR values of the S — R and R — D links, respectively. When m SR = 
mRD = m, the distribution in (|36l ) is reduced to 



(7717)™ _( 1 + 1 ) mi 



f fry) — - e VTSfl ' TRD 

/7 ™" l7J r(m) 7 



x 



1 1 / 7717 \ | 1 >p 1 / r»7 \ 

•R) m ~ kl \1rdJ (lR D ) m t-t k\ \jsr) 



(37) 



The average value of (T331) when m SR = m RD = m 



7 mm = -2(-l)- m ^l ( 7/?D B (-^-,m + 1, -2m 



r(m) 2 I V 7 flZ) 



+ 7sflB ( -^P^,m+ 1, -2m 

75ii 



(38) 



where B(x,y, z) is the incomplete Beta function. 

The minimum equivalent relay of the N-source relay channel in Fig. [2] has the instantaneous 
SNR 

7 min = min{7s li? , 7s 2jR , . . . ^ Sn r, Ird} (39) 

where jSnR is the instantaneous SNR per symbol of the channel S n — R, 1 < n < N. The 
distribution of (l39l) and its average SNR value, however, is not as readily found for large values of 
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N, making the model unsuitable for error rate analysis but useful for designing a low-complexity 
receiver, which is described below. 

The minimum equivalent hard-receiver of the relay network in Fig. [T] obtains the ML estimates 
of Ui by the hard-decision ML rule in (fl"5l) , where the SEP terms of the network coded data in 
the likelihood function in (1271) are approximated to the coherent M-PSK demodulation SEP of 
a minimum equivalent relay channel which has an instantaneous SNR value per symbol 7 m j„: 

P min = Q(V 2 lmin) + ~ / 6Xp - (it - ^^minf Q [V^U tan(7r/M ) J du 

«-i '" (40) 

= 2j V (zd,tic\ur = u R + e R ) p(e R ). 

An approximate error rate of the relay network in Fig. \T\ is found using an equivalent relay 
channel model in the following manner. The channel over which the network coded data at the 
terminals are sent to the destination are approximated using an equivalent channel model. The 
average PEPs of the union bound in (fTTT) are approximated using the same method employed in 
the error-free case, given by (fl"8l) . except that the SNR vector is obtained using the equivalent 
relay channels. 

Let us illustrate how the Q-inverse and minimum equivalent relay channels are used in 
obtaining the PEPs with an example. Consider the 2-user GF(4) NC, DetF, QPSK modulated, 
relay network described in previously in (|29l)-(l34l) assuming a Nakagami-m fading environment 
with m=l. The relay network has a generator matrix given by (|29| ) and the observed SNR 
vector at D is 

r = [7id ,sys l2D,sys 1eq,l 1eq,2\ 

(41) 

where r y n D,aya is the instantaneous SNR value per symbol of the channel over which u n is 
transmitted and 7 eg n is the instantaneous SNR value per symbol of the equivalent relay channel 
over which u nR is sent, n — 1,2. Assuming the S — R and R — D channels have equal average 
SNR values, the average values of fiD are ^i D , sy s = l2D, sys = 7 and = % gt2 = % q and 
the corresponding average SNR vector is 

f=[7 7 % q % q \ (42) 



Inserting m — 1 and ^sr — Ird = 7 in (1381) . the minimum equivalent relay channel's average 
SNR value is found as ^f m in —7/2, hence the average SNR vector when the minimum equivalent 
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relay channel model is used is the following: 

f min =[7 7 7/2 7/2]. (43) 

The average SNR vector when the Q-inverse equivalent relay channel model is used is given 
below: 

f q i = [7 7 7Q- 1 Iq-A ( 44 ) 

where 7Q 1 is found by inserting (1321) and (|34|) into ([27]) for M = 4. The average PEP between 

the codewords X(0) = <p([0 0]) and X(l) = <p([0 1 2 1]) is found using the MGF 



method in (fl~8l): 



1 , >(o) -^(o)| 2 \ _1 / , j^(o)-^(i)f ' 



P(X(0HX(1)) =Z / ( l + 7 - — 1 + 7 w 



7T 

6»=0 

-2 



Mi ^:5 i « (45) 



4sin 2 (#) 

1 7V 2sin 2 (£?) \ / 2sm 2 (#) 



n J \2 sm\6) + 7/ V2 sin 2 ((?) + 



d6 



where the modulated codewords are taken as <^(0) = 1, (p(l) = j, (p(2) = —j and 99(3) = — 1, 
assuming Gray coding. When the minimum equivalent relay SNR vector in (1431 is used in (|45l) , 
the PEP reduces to 

1 7V 2sin 2 (0) \ / 4sin 2 ( L 9) x 2 



P X(0W X(l) = - / 7 ; 7 I; ; d0. (46) 

V yJ y! ) nj V2sin 2 (#) +7; V4sm 2 (^) +7 ' 

6»=0 

Inserting the Q-inverse equivalent SNR vector in (|44|) into (l45l) gives 

p (x(oh x(D) - i /V 2 t? } 1 ( \f {B) ) dfl. (47) 

V KJ K> ) tt 7 V2sin 2 (#)+7A2sin 2 (#)+7 Q -iy' 

For high SNR values, using the fact that the Q-inverse equivalent relay channel converges to 
the minimum equivalent relay channel ifTTl . the union bound from (1461) and (1471) is found using 
Formula 3.621.3 on p.389 in IfTTl : 

tt/2 

p ( x( o H x (1) ) < -| / .AW = ~ = yr m 

8=0 

It is seen after this analysis that the diversity order in this case is equal to 3, from (|48l ). 
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IV. Performance Results 

In this section, optimal receiver's error performance of the relay network shown in Fig. Q] is 
compared to the error performances of the minimum and Q-inverse receivers and the approximate 
union bounds obtained using the equivalent relay channel models. For each SNR value, at least 
50 bit/symbol errors were observed in obtaining the corresponding BER/SER. 




LU 10 

m 1 



2 users, optimal hard-receiver 
2 users, Q-inverse eq. hard-receiver 

2 users, minimum eq. hard-receiver 

3 users, optimal hard-receiver 
3 users, Q-inverse eq. relay 

3 users, minimum eq. hard-receiver 



10 

SNR (dB) 



Fig. 5. 2 and 3-user GF(2) NC DetF Nakagami-m (m=l) fading relay networks' hard-decision ML performances 



Fig. [5] compares the optimum hard-receivers' BER performances of the 2- and 3-user cases 
of the GF(2) NC DetF Nakagami-m fading with m = 1 (Rayleigh flat-fading) relay network in 
Fig. Q] with that of the minimum and Q-inverse equivalent relay channel receivers. The generator 
matrices of the 2- and 3-user relay networks are given in (l49l) and (|50l) . respectively. 

10 11 
111 



G 2 



(49) 



G, 



I 1 1 1 1 \ 
10 110 
1 1 1 j 



(50) 
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Optimal hard-receiver, m=2 

-* — Minimum eq. hard-receiver, m=2 
-9 — Q-inverse eq. hard-receiver, m=2 
-<l — Optimal hard-receiver, m=1 

— Minimum eq. hard-receiver, m=1 
-a — Q-inverse eq. hard-receiver, m=1 



10 



SNR (dB) 



15 



Fig. 6. 2-user GF(2) NC DetF Nakagami-m fading relay network for m=l and m=2, hard-decision ML performance 
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Fig. 7. 2- vs 3-user GF(2) NC DetF Nakagami-m (m=l) fading relay networks' soft-decision ML performances 
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- 2 users, optimum soft-receiver 






- 2 users, minimum eq. relay union bound 




— e- 


- 2 users, Q-inverse eq. relay union bound 






- 3 users, optimum soft-receiver 




X — 


- 3 users, minimum eq. relay union bound 




— e— 


3 users, Q-inverse eq. relay union bound 
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Fig. 8. 2- and 3-user GF(2) NC DetF Nakagami-m (m=l) fading relay networks' optimum soft-receiver BER performances vs 
minimum and Q-inverse equivalent relay networks' union bound BER performances 



The Q-inverse and minimum equivalent relay channel hard-receivers are observed to display 
an error performance very close to each other and to the optimal hard-receiver. The diversity 
order is shown to increase with an increasing number of users; the diversity orders of the 2 and 
3-user relay networks are 2 and 3, respectively. Fig. |6] compares the BER performance of the 
optimal hard-receiver of the 2-user GF(2) NC DetF Nakagami-m fading relay network in Fig. \T\ 
with that of the minimum and Q-inverse equivalent relay channel hard-receivers for m = 1 and 
m = 2. The generator matrix is given in (|49l ). It is observed that when m = 2, the minimum and 
Q-inverse equivalent channel receivers display a very close error performances to the optimal 
receiver and have a diversity order close to 3, compared to a diversity order of 2 when m = 1. 

Fig. U} compares the optimum soft-receiver BER performance of the of the 2- and 3-user GF(2) 
NC DetF Nakagami-m fading relay network in Fig. \T\ with that of the minimum and Q-inverse 
soft-receivers when m — 1. The generator matrix used is given in (|49l) and (|50l) . respectively. 
It is observed that the Q-inverse equivalent relay channel soft-receiver displays a very close 
error performance to the optimal soft-receiver, whereas the minimum equivalent relay channel 
soft-receiver requires 1.7 dB less SNR for a BER=1CT 5 in the 2-user case and a diversity order 
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Fig. 9. 2-user GF(4) NC DetF Nakagami-m (m=l) fading relay network's optimal soft-receiver SER vs minimum and Q-inverse 
equivalent relay networks union bound SER performances 

loss in the 3-user case. 

Fig. [8] compares the analytical union bound BER performances of the 2- and 3-user GF(2) 
NC DetF minimum equivalent and Q-inverse equivalent relay networks with the optimum soft- 
receiver BER performances of the 2- and 3-user relay networks in a Nakagami-m fading en- 
vironment when m = 1. The relay networks' generator matrices are given in (|49|) and (|50l) . 
respectively. The minimum and Q-inverse equivalent relay networks' union bounds are very 
close to the relay network's optimal performance. 

Fig. M compares the analytical union bound SER performances of the 2-user GF(4) NC DetF 
minimum equivalent and Q-inverse equivalent relay networks with the optimum soft-receiver's 
SER performance of the relay network in a Nakagami-m fading fading environment when m — 1. 
The relay network's generator matrix is given in (|29l) . The union bound SER performances of 
the Q-inverse and minimum equivalent relay networks are very close to each other and approach 
the optimum soft-receiver SER performance at high SNR values. 

Fig. [10] compares the analytical union bound BER performance of the 2-user GF(2) NC DetF 
Q-inverse equivalent relay network with the optimum soft-receiver's BER performance of the 
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Fig. 10. 2-user GF(2) NC DetF Nakagami-m (m=2) fading relay network's optimal soft-receiver BER vs Q-inverse equivalent 
relay network's union bound BER performance 

relay network in a Nakagami-m fading fading environment when m — 2 . The relay network's 
generator matrix is given in (l49i The union bound BER performance of the Q-inverse equivalent 
relay network is very close to the optimum soft-receiver BER performance. 

V. Conclusion 

This paper presents simple means of analyzing the error rate performance of a general wireless 
cooperative GF(q) NC DetF relay network with known relay error statistics at the destination 
using equivalent channel models. The approximate error rate of the relay network is derived with 
the use of equivalent relay channel models, where the minimum and Q-inverse equivalent relay 
channel models are adapted to N-user GF(q) NC DetF Nakagami-m fading relay channels. The 
equivalent channel models are used in developing high-performance ML receivers for Nakagami- 
m fading channels, which are shown to display very close error performances to the relay 
network's optimum receiver. 
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